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Abstract 

Introduction. One of the reasons for undesirable delamination of polymer composites with fabric reinforcement is low 
transverse shear properties. It is known that the reinforcement of polymer fabric composites in the Z direction reduces the 
sensitivity to delamination and increases the viscosity of interlayer fracture. Various methods of three-dimensional 
reinforcement of polymer fabric composites are proposed in the literature. However, they complicate the manufacturing process 
of the structure. The problem is solved by the method of three-dimensional reinforcement proposed in this article — felting. 


This is a local reinforcement of the composite in the Z direction with minimal production changes. The degree of 
Z-reinforcement is determined by the felting density, i.e., the number of needle punches per 1 cm? of the fabric package. The 
work is aimed at evaluating the effect of felting on the interlayer crack resistance of a composite material. 

Materials and Methods. The interlayer fracture toughness Gye was determined on a cross-woven fiberglass with 
felting of 10 cm?. The material was impregnated with Etal-370 resin and Etal-45 hardener. Experiments according to 
ASTM D7905M-14 and GOST 33685-2015 standards were carried out on an Instron 5900R test machine. The stress 
state at the crack tip was analyzed with regard to the nonlocal strength theory in the ANSYS Workbench program (option 
“static strength analysis”). The finite element method (FEM) was used. 

Results. The “load — displacement” curves were considered for the samples. Values Gy. were calculated. The results of 
EMF tests for felting density of 0 cm and 10 cm were summarized. Control samples and felting samples were compared. 
In the latter case, Gu. turned out to be ~33% higher. The stress state at the crack tip was calculated under DCB and ENF 
loading. The dependences of maximum normal and shear stresses, as well as displacements, were visualized in the form 
of graphs and color charts. To get the calculated “load — displacement” dependences using FEM, the reverse method of 
obtaining transverse shear constants was used. DCB loading showed that felting provided increasing the rupture strength 
in the Z direction to ~18%, by 39 to 46 MPa, and in the planes XZ— to ~16%, by 77 to 89 MPa. 

Discussion and Conclusion. Felting as a method of local three-dimensional reinforcement enhances the interlayer crack 
resistance of polymer fabric composites. It provides reducing the area of stratifications after local impacts during the 
operation of structures. Flexible felting technology makes it possible to create zones with an arbitrary impact density, 
increasing fracture toughness only in the required places of structures. The FEM analysis of the stress state at the crack 
tip within the framework of the nonlocal strength theory has shown that in strength calculations, the stratification crack 
can be considered as a stress concentrator. 


Keywords: reinforcement of polymer fabric composites, transverse shear strength, interlayer crack resistance, 
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AHHOTalna 


Beedenue. O7Ha 43 NpW4WH HeKeNAaTebHbIX PaccOCHHM NOJMMepHbIX KOMIMO3HTOB C TKAHEBbIM apPMHpOBaHueM — HH3- 
KWe TpaHCBepCasIbHO-CABHTOBbIe XapakTepHcTHKH. M3BecTHO, YTO aPMMpOBaHHe MOJIMMepHbIX TKAHEBbIX KOMIO3HTOB B 
HarpaBuleHHu Z yMeHBIULaeT YYBCTBUTEJIBHOCTh K pacCJIOCHHI0 H MOBbILMAeT BASKOCTbh MEXKCIIOHHOTO paspymieHua. B smTe- 
paType lpeyyiararorca pa3Hble ClocoObl TpexMepHOrO aPMHPOBaHHA MOIMMePHbIX TKAHEBbIX KOMIIO3HTOB. OZHaKO OHM 
YCJIOXKHAIOT IIpouecc H3rOTOBIeHHA KOHcTpyKuMH. IIpoOmemy peliaeT IpesIO%KeHHEIM B JaHHOM CTaTbe clOcoO TpexMep- 
HOrO apMupoBaHna — dbeuITHHT. ITO JOKabHOe apMHpOBaHHe KOMMO3HTa B HallpaBIeHHH Z IIpH MHHMMAaJIbHBIX MIpOH3- 
BOJICTBCHHBIX H3MeHeHHAX. CTerleHb Z-apMupoOBaHHA OMpeeAeTCA TWIOTHOCTHIO (pesITHHTA, T.e. KOJIMYCCTBOM y{apoB 
urbI Ha | cM? rKanesoro maxkerta. [ens paSoTbl — OLeHHTh BIIMAHHe (eITHHTAa Ha MEXKCIIOMHY!O TpelIMHOCTOMKOCTS KOM- 
MO3HTHOrO MaTepnalia. 

Mamepuajei u memoooi. MexcnoiHy0 BaA3KOCTb pa3pyleHHa Gy. OlIpeeMAIM Ha CTCKIIOTKAHM MOJIOTHAHOYO Mepe4e- 
TeHua c (pertunrom 10 cm’. Marepuam mponmtTsisaru cmMomoi Sran-370 u orBepaqurenem Sran-45. SkcrepuMeHrTs! 110 
cranfaptam ASTM D7905M-—14 u OCT 33685-2015 nposoyumm Ha HcnbItTaTesIbHOM Mamuue Instron 5900R. Hanpsxen- 
HOe COCTOAHHE Y BEPLIMHbI TPeLMHbI aHaIM3HPOBAIM C MO3HI[MM HeOKaJIbHOM TeOpHu NpOuHOCTH B IporpamMe Ansys 
Workbench (onmma «cTaTH4ecKHi MpOUHOCTHOM aHas3»>). 3a\eHCTBOBAaJIM MeTOJ, KOHCYHBIX 9IEMeHTOB (MK3). 
Pezynomamoi ucciedosanua. [Ina oOpa3yoB paccMOTpeM KpuBble «Harpy3Ka — TepeMeljeHHe». Baranciuiu 3Ha4e- 
Hua Gye. O60OujuIu MToru ENF-ncnpiranuii W14 WOTHOCTH (berTuHTa 0 cm 4 10 cm. CpaBHusIM KOHTpOUbHEIe OOpasiib 
WM OOpa3ubI c (berTHHToM. B mocueyqHeM ciy4ae Gre OKa3aacb BbILIe Ha ~33 %. PaccuuTaM HallpsKeHHOe COCTOAHHE y 
BePLIMHBI TpewHHb! mpx DCB- u ENF-varpyxennu. Busyarm3npoBasin B BAe rpa:bHKOB H I{BeTOBbIX {MarpaMM 3aBHCH- 
MOCTH MaKCHMA€JIbHBIX HOPMA@JIBHBIX M KaCaTeJIbHBIX HallpsKeHHH, a TakKoKe NepememjeHni. JI nowy4aeHuad pacueTHBIX 3a- 
BHCHMOCTeli «Harpy3Ka — TepemeljeHve» c TomMonbIo MKD ucnomb30BaiM OOpaTHbIM MeTOA MouyYeHuA TpaHcBep- 
CaJIbHO-CA/BHTOBBIX KOHCTaHT. Harpy2keHue 110 cxeme DCB noxa3aui0, 4TO (PeITHHT MO3BOJACT YBEIMYHT pees IpOuHo- 
CTH Ha pacTsKeHue B HallpaBsieHun Z Ha ~18 %, c 39 Wo 46 MIa, a B mIocKoctu XZ — Ha ~16 %, c 77 MI a no 89 MIa. 
O@cystcoenue u 3aknio4“enue. DeNTHHT Kak CIOCOO JOKaIbHOrO TPeXMePHOrO apMUPOBaHHA YCHIMBaeT Me2KCIIOMHYIO 
TPeLWIHHOCTOMHOCTh MOJIMMepHbIX TKAHEBbIX KOMIO3HTOB. OH NO3BOJIAeT COKPATHTb WIOWayb paccioeHu Mocse J0- 
KaJIbBHbIX Y1apoOB IIpH SKCIUIyaTaljuu KOHCTpyKuM. [uOKad TeXHONOrMA (PeNTHHTAa TaeT BOSMOXKHOCTL CO3LaBaTb 30HbI 
C MIpOH3BOJbHOM MWJIOTHOCTEIO YapoB, MOBbIMMAA TPeLIMHOCTOMKOCTb JIMIIb B HeOOXOAMMbIX MecTaxX KOHCTpyKUHH. 
MK®9-aHasIM3 HallpsoKeHHOTO COCTOAHHA Y BEPLIMHbI Tpel{HHbI B paMKaX HEOKAIbHOM TEOPHU NPOUHOCTH NOKa3all, YTO 
B IIPOUHOCTHBIX pacueTax TpelHHY PaccIOCHHA MOXKHO paCCMaTPUBATb KaK KOHICHTpaTOp HallpsKeHHH. 


Kunrouesble CJ10Ba: apMUpoBaHHe MOJIMMeCpHBIX TKAHCBEIX KOMHO3HTOB, TpaHCBepCaJIBHO-CABHTOBaxA TIPOUHOCTh, MexKcowHad 
TPeLIMHOCTOMKOCTS, MO?*KCJIOCBAA BAZSKOCTb pa3pyuieHHaA, (pesITHHTOBOe JIOKaJIBHOe TpexXMepHoe apMHpoBaHHe 
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Introduction. Fibrous polymer composites are widely used, in particular, in aviation and space engineering due to 
their significant rigidity and strength in the fiber orientation (plane XY) [1]. However, the transverse shear strength of 
these materials is quite low [2], as it is determined by the features of the polymer matrix [3]. Reinforcement of polymer 
fabric composites in the Z direction provides reducing the sensitivity to delamination, i.e., increasing the viscosity of 
interlayer fracture [4]. 

Various methods of three-dimensional reinforcement of polymer fabric composites are known [5]. However, they 
create additional difficulties in the manufacture of structures made of polymer fabric composites [6]. The method of three- 
dimensional reinforcement proposed in this work, felting [7], makes it possible to obtain a locally reinforced composite 
in the Z direction with minimal changes in the production process. The degree of Z-reinforcement is determined by the 
felting density, i.e., the number of needle punches per | cm? of the area of the fabric package [8]. 

The use of various methods for determining the fracture toughness of polymer composite materials [9] makes it 
possible to conduct studies on various samples [10] and with different loading methods [11]. One of the most common 
approaches is the three-point bending method. In this case, a beam-shaped delamination sample is used. We are talking 
about ENF tests (End-Notched Flexure — bending of a sample with edge delamination) [12], which involve transverse 
shear loading. This makes it possible to determine the interlayer fracture viscosity of Ge — mode II fracture. High shear 
stresses occur at the crack tip [13]. 

Another common way to determine transverse characteristics is the Double Cantilever Beam method (DCB tests). In 
DCB tests, the value of the interlaminar fracture toughness Gj). is determined under separation loading — fracture 
according to mode I [14]. The delamination crack spreads due to the action of normal stresses [15]. 

The presented work was aimed at the evaluation of the effect of felting on the interlayer fracture toughness of a 
composite material. To do this, ENF tests (bending of a sample with an edge separation) of a composite material with 
increased crack resistance due to the use of felting were carried out. Previously, the authors studied the effect of felting 
on the interlayer crack resistance of a composite material during DCB tests [16]. A computational model based on the 
nonlocal theory of strength has been developed. It provides for the calculation of the stresses that occur in ENF and DCB 
samples, for cracks of various lengths, using the finite element method (FEM). 

Materials and Methods 

Experimental determination of interlayer crack resistance by the ENF method. The samples were made of cross- 
woven fiberglass with a layer thickness of 0.2 mm. A package of dry two-layered fiberglass was punched on a felting 
machine with a felting density of 10 cm (10 punches with a felting machine needle per 1 cm? of dry glass fabric). The 
fiberglass package was punched in such a way that after impregnation and hardening, the initial crack did not fall on the 
felting area. To create an initial crack between two layers of fiberglass, an aluminum foil with a thickness of 11 um, 
coated with a Vs-M parting lubricant, was placed. Fiberglass was impregnated with resin Etal-370 and hardener Etal-45. 
For the production of reference samples, two layers of dry fiberglass were impregnated with Etal-370 resin and Etal-45 
hardener without punching on a felting machine. After impregnation, plates from glass textolite STEF (electrotechnical 
glass-cloth-base laminate) were glued to two layers of fiberglass (Fig. 1 a). Samples with a length of 150 mm and a width 
of 16 mm were obtained by cutting the hardened plates using a high-speed circular saw. 


Fiberglass STEF 2 mm 


Felting 
area 


a) 


Fig. 1. Configuration and parameters of ENF tests according to mode J/: 
a — loading scheme; b — photo of tests 
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The Instron 5900R test machine with a loading speed of 10 mm/min was used. The distance between the supports was 
2L = 100 mm. The initial crack length for all samples was ap = 25 mm. To perform compliance calibration over a wide 
range of crack lengths, one sample of each type was unloaded and reloaded. The resulting crack served as the initial crack 
in the next loading cycle. 

Calibration recommended by ASTM D7905M-14! and GOST 33685-2015? standards was used to process the test 
results. This approach provided determining parameters A and m for each felting sample and each non-felting control 
sample from the linear dependence of two quantities — compliance of sample C and cube of the crack length a?: 


C(8/P(8))=A+m-a’, (1) 
where P — load applied to the sample; 6 — displacement. 
After calibration and determination of parameters A and m, the crack length can be found from expression (1): 


o-(o=4)". re 


m 


The moment of the delamination onset is determined by the condition C(6) = C(Pmmax). Value of the interlayer fracture 
toughness at the separation onset (crack development): 


3m: Po a? 
Gite = —— (3) 
where Pmax — maximum load; a — crack length calculated by formula (2) at load Pina; b — sample width. 

Calculation of the stress state at the crack tip under loading according to the DCB and ENF schemes. The stress 
state of a crack-like concentrator is estimated from the perspective of approaches that use nonlocal stresses [17], i.e., 
averaged on some basis [18]. The calculation model also includes the assumption of linear-elastic behavior of the material 
up to destruction [16]. 

The main hypothesis is that the strength criterion of the composite, which includes all components of stress averaged 


on the basis A, is responsible for the development of a crack-like concentrator (Fig. 2): 


2 2 2 
max O Max O >» Max T MAX O ,, * Max oO 
2] { 2) +( =| zh XA <1, (4) 


Z, Xx, Ss Z,°X;, 
where Z, and _X;— rupture strength in the Z and_X directions; S — shear strength in the plane XZ. 


— 
{ Max On 4 MAX Ox, 
a) b) 


Fig. 2. Stresses averaged on the basis A at the crack tip: 
a— ENF tests; b — DCB tests 


' ASTM D7905/D7905M-14. Standard Test Method for Determination of the Mode II Interlaminar Fracture Toughness of Unidirectional Fiber-Rein- 
forced Polymer Matrix Composites. URL: https://cdn.standards.iteh.ai/samples/89096/03be6b5e53664f13a8703bb4342d98 la/ASTM-D7905- 
D7905M-14.pdf (accessed: 22.04.2024). 

2 GOST 33685-2015. Polymer Composites. Test Method for Determination of the Interlaminar Fracture Toughness under Shear. (In Russ.) URL: 
https://docs.cntd.ru/document/1200127774 (accessed: 22.04.2024). 
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Due to the presence of symmetry planes, a three-dimensional 1/2 crack sample model was constructed for ENF 
loading (Fig. 3), and 1/4 crack sample — for DCB tests (Fig. 4). Calculations were performed in the ANSYS Workbench 
program (option “static structural’’). 


Contact 
Frictionless 


Contact 
Bonded 


Fig. 3. Finite element model 1/2 of the sample and a fragment of the grid for ENF tests 


0.0 30.0 60.0 (mm) 


= SC 


15.0 45.0 


Fig. 4. Finite element model 1/4 of the sample and a fragment of the grid for DCB tests 


When creating the finite element grid, parameter 1, = 0.75 mm in width of the sample [19] and parameter 1, = 0.2 mm 
in thickness of the sample were set, which corresponded to the thickness of the modified layer [20]. One finite element in 
the thickness of the layer was set in accordance with the layer wise theory used in assessing the strength of layers within 
the framework of mesomechanics of composites [21]. In ENF tests, the total displacements in the sample were much 
greater than the local displacements from the loading roller (Fig. 1 5); therefore, the grid of finite elements was not 
condensed at the places of application of loads and supports (Fig. 3). Properties of fiberglass used in the calculation: 

- elastic modules E,. = Ey = 23 GPa, E, = 9 GPa; 

- shear modules Gy, = Gy: = Gy = 6000 GPa; 

- Poisson's coefficients [x = 0.15, Uy: = Hx: = 0.3 [22]. 

Since the volume fraction of transverse reinforcement is less than 1% [16], it is assumed in the calculations that the 
elastic properties of fiberglass do not change under felting. 

Dependence P(6) was calculated in accordance with the sequence described below. 

1. FEM-calculation of the maximum stresses maxoy, maxoy and maxtx, and displacements of point 6 of application 
of load P= 1H for cracks with the given lengths in the range a = 20...90 mm (DCB) and a = 25...40 mm (ENF) was 
performed. 

1. Approximation dependences o. = fla, P) = P:by-a; 6, = fla, P) = P-bra; %a.= fla, P) = P-b3:a; 5 = fla, P) = P-cy-a*-(DCB) 
and oa = fla, P) = P-by-a; 6n.= fla, P) = P-bra; twa = fla, P) = P-(bs-a + ds); 8 = fla, P) = P(cv-a? + cx.a® + c3:a + c4)(ENF). 

2. were constructed using the least squares method. 


3. Load P.,(ao) and displacement 6,,, at which the initial crack length ao would increase abruptly by Ax = 0.75 mm 
upon the violation of the strength criterion, were determined (4). 
4. With crack length agtnA, loads P(aotnA) and displacements 6(n) for n>0 were determined. 
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Research Results 

Results of the experimental determination of interlayer fracture toughness by the ENF method. Figure 5 shows 
the “load — displacement” curves for all tested samples. All “load — displacement” curves have an area with constant 
compliance (C;;,), corresponding to the linear “load — displacement” ratio. Values Cj, were used for calibration. 


P,kKN 


Fig. 5. “Load — displacement” diagrams of ENF tests 


without felting; = = = =— with felting 


Figure 6 shows the calibration curves. For felting and non-felting samples, compliance is proportional to the cube of 
the crack length. 


C, mm/N 
0.012 ---O 
y=8.44- 10°. + 6.69: 10° el _*2 
0.010 ae pt ee 
pot a 
0.008 =  -8- 
goer? y=6.98-10°x+6.29-10° 
0.006 
0.004 
0.002 
0.000 
0 10,000 20,000 30,000 40,000 50,000 60,000 a, mm? 


Fig. 6. Dependence of compliance of sample C on the cube of crack length a?: 
o — without felting; e — with felting 


To calculate values of the crack length a*, corresponding to the compliance at the delamination onset C(Pinax), the 
obtained calibration curves and equation (2) were used. For the found values of the crack length a*, values Gy. were 
calculated using equation (3). The results are shown in Table 1. 
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Table 1 
ENF Test Results 

ee ao, mm | Cin, mm/N ear a*,mm | Pax, N Gite, kI/m? ae ao 
0 25 8.032 8.845 29.4 504.6 1.723 
0 25 7.980 8.807 29.3 544.6 1.983 
0 25 8.299 9.392 31.7 489.4 1.908 pe One 
0 25 7.905 8.900 29.7 505.2 1.746 
10 25 7.587 8.337 30.8 627.3 2.432 
10 25 7.849 8.677 32.4 625.8 2.682 
10 25 7.937 8.811 33.0 578.4 2.376 2.441 + 0.154 (6.3%) 
10 25 7.824 8.594 32.1 581.4 2.261 
10 25 7.880 8.818 33.1 589.3 2.456 

* variation coefficient 


Felting samples showed a significant (by ~33%) increase in the interlayer fracture toughness Gy. compared to the 
control samples. After testing, felting samples were separated with a sharp knife and examined under a microscope. 
Micrographs of the zone without felting (area of the initial crack) and the zone with felting (area of crack development) 
are shown in Figure 7. When cracks develop, the fibers elongated during felting are destroyed, because their length is 
greater than the critical one [16]. 


Fig. 7. Micrographs of felting samples after ENF tests: 
a — zone without felting (area of initial crack); 


Mechanics 


b — zone with felting (area of crack development); 
c — zone with felting (enlarged scale) 
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Results of the calculation of the stress state at the crack tip under loading according to the DCB and ENF 
schemes. Figures 8—9 show the dependences of stress maxoa(a), maxow(a), maxtxx(a) and displacements 8(a). 
Conditions: P = 1 N, loading according to DCB and ENF schemes. 


max 6,, MPa max 6,, MPa 1 
6.0 y=2.07- 10x 
1.6 
5.0 R’=0.99 
4.0 1.2 
3.0 08 
2.0 
1.0 0.4 
0.0 0.0 
0 20 40 60 80 a,mm 0 20 40 60 80 a,mm 
a) b) 
MAX Ta, MPa 6,mm 
0.6 1.2 
0.4 0.8 
0.2 0.4 
0.0 0.0 
0 20 40 60 80 a,mm 0 20 40 60 80 a,mm 
c) d) 


Fig. 8. DCB-loading. Dependences of maximum stresses and displacements of the crack length at P= 1 N: 
a— dependence of normal stresses maxo.a(a); b — dependence of normal stresses maxoa(a); 
c — dependence of shear stresses maxtxa(a); d — dependence of displacements 5(a) 


max 6,, MPa >, —_ —___,_—__"——_ max 6;, MPa 
y=4.13-10°x y=4.13- 10x 
0.20 0.20 
R’ =0.999 R’ =0.999 
0.15 0.15 


0.00 + 0.00 + 
0 10 20 30 40 a,mm 0 10 20 30 40 a,mm 
a) b) 
Max T%,, MPa A 6, mm 1 ch, 5 : r 
y=-9.40-10 x + 1.63 -10°x — 4.63 - 10° x+ 
0.20 
+ 1.09- 10° 
0.15 0.12 ; } 
2 R= ae 
8 0.10 0.08 | —, 
=| 
3 0.05 0.04 | 
pI 
& 0.00 4 0.00 
g 0 10 20 30 40 a,mm 0 20 40 60 80 a,mm 
— 
Pa c) d) 
b= Fig. 9. ENF- loading. Dependences of maximum stresses and displacements of the crack length at P= 1 N: 


a — dependence of normal stresses maxo.(a); b — dependence of normal stresses maxoz(a); 
c — dependence of shear stresses maxtxa(a); d — dependence of displacements 8(a) 
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Examples of stress distribution at the crack tip are shown in Figures 10—11 with crack length ap = 30 mm. 
aS a 


2.108 1.325 0.361 0.512 —1.822 MPa 0.639 0.450 0.262 0.073 -0.115 MPa 
a) b) 


aS Eee 
0.364 0.234 0.104 -0.026 -0.156 MPa 


¢) 


Fig. 10. Stresses at the crack tip under DCB loading: a — normal stresses oy; 
b — normal stresses oz; c — shear stresses Txza 


0.649 0.361 0.072 —0.216 -0.505 MPa 0.053 0.292 0.006 —0.018 —0.041 MPa 


a) b) 
a | 


0.155 0.112 0.069 —0.026 —0.017 MPa 
°) 


Fig. 11. Stresses at the crack tip under ENF loading: a — normal stresses ox; 
b — normal stresses 6; c — shear stresses Tx 


To obtain the calculated “load — displacement” dependences using FEM, the strength characteristics of the composite 
in the main directions are taken into account, i.e., criterion (4). Direct obtaining of transverse shear constants involves a 
certain difficulty; therefore, the reverse method is considered below. With this approach, the constants vary, and their 
best combination is found. This means that the calculated and experimental loading diagrams are in good agreement 
(mean square deviation of displacements at specified loads is minimal). 

The results of the calculation under loading according to the DCB scheme were compared to the authors' experiment, 
which was considered in [16]. Samples were made in the same way. The tests were carried out in accordance with GOST 
R 56815-20153 and ASTMD5528-14¢ standards. 

The calculation was performed for loading according to the DCB scheme. When calculating dependence P(S) for 
samples without felting, the following stress limits were found and rounded to integer values: Z;= 39 MPa, X; = 360 MPa, 
S = 82 MPa. The obtained values 360 MPa and 39 MPa correspond to the data on the strength of fiberglass specified 
in [23]. For felting samples (density 10 cm”), the calculated values were X;* = 270 MPa, Z;* =46 MPa and S* = 97 MPa. 
Thus, the use of felting made it possible to increase the rupture strength in the Z direction from 39 to 46 MPa (by ~18%). 


3 GOST R 56815-2015. Polymer Composites. Method for Determination Specific Work of Exfoliation in Tearing Off Conditions. (In Russ.) URL: 


https://docs.cntd.ru/document/1200131393/titles (accessed: 22.04.2024). 
4 ASTMD5528M-21. Standard Test Method for Mode I Interlaminar Fracture Toughness of Unidirectional Fiber-Reinforced Polymer Matrix Compo- 


sites. https://doi.org/10.1520/D5528 D5528M-21 
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When loading according to the DCB scheme, the shear strength limits in the plane XZ S and S* do not make a big 
contribution to criterion (4); therefore, values S = 82 MPa and S* =97 MPa obtained in calculations according to the 
DCB scheme need to be clarified according to the ENF loading scheme. Note that the effect of normal stresses in the X 
and Z directions is insignificant compared to shear stresses under loading according to the ENF scheme. Therefore, in the 
calculations, when searching for the values S and S*, values Z,, X;, Z;* and_X;* were taken from solving the inverse problem 
under loading according to the DCB scheme. 

Values S = 77 MPa (without felting) and S* = 89 MPa (with felting) were determined from the condition of the best 
consistency of the experimental and calculated curves P(5) (mean square deviation of displacements at given loads is 
minimal). Evidently, felting made it possible to increase the shear strength in the plane XZ by ~16%. 

Figures 12-13 show the experimental “load — displacement” diagrams, as well as calculated dependences P(6) for 
the found values of stress limits: 

- Z,;= 39 MPa, X; = 360 MPa, S= 77 MPa (for control samples without felting); 

— Z;* = 46 MPa, X;* = 270 MPa, S* = 89 MPa (for samples with felting). 


Fig. 12. Experimental diagrams “load — displacement” [16] and calculated dependences P(8): 

— DCB tests of samples without felting; — calculation without felting (Z;= 39 MPa, X;= 360 MPa, S = 77 MPa); 
- -- DCB tests of samples with felting; - - - calculation with felting (Z;* = 46 MPa, X;* = 270 MPa, S* = 89 MPa) 
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Fig. 13. Experimental diagrams “load — displacement” and calculated dependences P(5): == ENF tests of samples without felting; 
— calculation without felting (Z;= 39 MPa, X;= 360 MPa, S= 77 MPa); - - - ENF tests of samples with felting; 
--- calculation with felting (Z;* = 46 MPa, X;* = 270 MPa, S* = 89 MPa) 


Discussion and Conclusion. Studies of the fabric composite have shown that felting with a density of 10 cm~ 
increases the viscosity of the interlayer fracture Gye by ~33%. 

Using FEM, the stress state was analyzed in a quasi-static elastic formulation of the problem and with a nonlocal strength 
theory for the developed numerical models of a beam with cracks of the known length. The distinctive feature of the 
calculations was that they did not use contact algorithms, but only considered the destruction of the composite layer closest 
to the crack, and the corresponding change in the area of gluing the layers. That is, the crack was considered as a stress 
concentrator. The composite strength criterion, which contained three parameters and was recorded through averaged 
stresses, provided using the method of step-by-step crack advancement to predict the “load — displacement” curve. 

The use of felting with a density of 10 cm™ increases the stress limit of the composite in the Z direction by ~18%, and 
the shear strength in the plane XZ — by ~16%. This became known from solving the inverse problem, i.e., searching for 
the strength characteristics of the material according to criterion (4) and the “load — displacement” curve. 

The results of the presented research will find their practical application. They can be specifically used in problems 
of forecasting defects, such as delamination (e.g., in low-speed impacts on composites in aircraft skin). The research 
results will be useful for eliminating these defects with the help of felting. 
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